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The so-called 'LSND anomaly', a 3.8a excess of Po events interpreted as originating from — > Pc oscillation, gave rise 
to many theoretical speculations. The MiniBooNE Collaboration reported inconsistency of this interpretation with 
the findings from their search for i/^ — > Ve oscillations. Yet the origin of the LSND anomaly was never clarified. A 
critical issue is the prediction of the background Pe flux that was used in the analysis of the LSND experiment. For 
this, decisive input comes from pion spectra measured with the HARP large-angle spectrometer under conditions that 
closely resemble the LSND situation: a proton beam with 800 MeV kinetic energy hitting a water target. 
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1. INTRODUCTION 
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The LSND experiment at Los Alamos studied the flux originating from protons with 800 MeV kinetic energy 
hitting a beam dump consisting of water to a good fraction. They reported a 3.8a excess of Dc events over back- 
ground The excess was interpreted as originating from — > oscillation, which gave rise to many theoretical 
I— i. speculations as to the existence of sterile neutrinos. 

^ ' The MiniBooNE Collaboration reported inconsistency of this interpretation with the findings from their search at 
I FNAL for oscillations 0. 

A critical issue in the LSND analysis is the background level of v^. events that originates from the decay chain 
,; 7r~ — > fjT — > Vg. An underestimate of 7r~ production which was quite uncertain at the time of the LSND experiment, 
would reduce the anomalous excess of Ve events. The relevant input is a precise measurement of the ratio of secondary 
TT^ to 7r+ production by 800 MeV protons. 

With a view to clarifying this issue, the HARP detector at the CERN PS that took data in 2001 and 2002 with 
0^ , proton and pion beams with momentum from 1.5 to 15 GeV/c, also recorded data from the exposure of a water 
target to protons with 800 MeV kinetic energy (1.5 GeV/ c momentum). 

This paper reports ratios 7r^/7r+ from this exposure and compares them with the ratios that were used in the 
LSND analysis. 
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2. DETECTOR CHARACTERISTICS AND PERFORMANCE 



The HARP detector combined a forward spectrometer with a large-angle spectrometer. The latter comprised a 
5^ ! cylindrical Time Projection Chamber (TPC) around the target and an array of Resistive Plate Chambers (RPCs) 
that surrounded the TPC. The purpose of the TPC was track reconstruction and particle identification by dE /dx. 
The purpose of the RPCs was to complement the particle identification by time of flight. For the work reported here, 
only the HARP large-angle spectrometer was used 0, 0| • Its salient technical characteristics are stated in Table HI 
The good particle identification capability stemming from dE/Ax in the TPC and from time of flight in the RPCs 
is demonstrated in Fig. [TJ Correct particle identification is of crucial importance as 800 MeV protons produce many 
more secondary protons than tt+'s. This is highlighted in Fig. [2] which demonstrates that the decomposition of the 
observed secondary particle spectrum into particle species is well understood. 
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Table I: Technical characteristics of the HARP large-angle spectrometer 



TPC 


RPCs 


cr(l/pT) ~ 0.20 - 0.25 (GeV/c)-^ 
<j{0) ~ 9 mrad 
a{dE/dx)/dEdxr.0.l6 


Intrinsic efficiency ~ 98% 
cr(TOF) ~ 175 ps 
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Figure 1: Specific ionization dE/dx (left panel) and velocity /3 (right panel) versus the charge-signed momentum of positive 
and negative tracks in -1-8.9 GeV/c data. 




Figure 2: Decomposition, on a logarithmic scale, of the observed spectrum of secondaries from the interaction of 800 MeV 
protons in water into particle species, as a function of the charge-signed 1/pt. 
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3. HARP-CDP MEASUREMENTS VERSUS LSND MONTE CARLO SIMULATIONS 

Most positive pions from the interaction of the 800 MeV protons in the beam dump are slowed down and decay 
at rest, while a few per cent decay at flight. The positive muons from the decay also slow down and mostly decay 
at rest, with only a very small fraction decaying in flight. The neutrino flux resulting from positive pions consists of 
z/^, Vfi, and v^. For negative pions the decay chain is the same except for charge conjugation, yet negative pions that 
come to rest disappear 100% by strong interaction, and likewise negative muons at the 90% level by weak interaction. 
The neutrino flux resulting from negative pions can only come from pions decaying in flight, and muons decaying 
in orbit after capture. The neutrino flux consists of D^^ Pc, and v^. Overall, the which is of interest is reduced 
by a factor of order 10"'*. In more detail, the reduction depends on the level of tt" production, the probability of 
7r~ decay in flight (related to the momentum spectrum and the beam dump geometry), and of the probability of a 
stopped /i~ to decay in orbit. 

Our measurement addresses the first of these three issues. Figure [3] shows preliminary results of the measured 
ratio of tt" to 7r+ production in four bins of polar angle 9 with respect to the incoming proton direction. Also shown 
is the parametrization of this ratio that was used in the LSND analysis Q . The measured ratio is smaller than the 
parametrization. Thus, if by the time of the LSND analysis our data had been known, the reported anomalous excess 
of Vc events would have been even larger. Further work on this interesting result is in progress. 




Pioo momBnlum (GeV/o) pion momentum (G«V/e) 

Figure 3: Ratios ■k~ /■k'^ in four bins of polar angle 6 with respect to the incoming proton direction, as a function of pion 
momentum; the crosses denote data, the lines represent the parametrization that was used by LSND. 
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